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Abstract Alfven waves are considered to be viable transporters of the non-thermal energy 
required to heat the Sun's quiescent atmosphere. An abundance of recent observations, from 
state-of-the-art facilities, have reported the existence of Alfven waves in a range of chro- 
mospheric and coronal structures. Here, we review the progress made in disentangling the 
characteristics of transverse kink and torsional linear magnetohydrodynamic (MHD) waves. 
We outline the simple, yet powerful theory describing their basic properties in (non-)uniform 
magnetic structures, which closely resemble the building blocks of the real solar atmosphere. 
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1 Introduction 

The solar coronal plasma is an ionized gas which is mainly confined by the ubiquitous 
presence of magnetic flux tubes and open magnetic field lines. These structures demonstrate 
sizes over very large scales, right down to the current observational limit, and are maintained 
at temperatures of several million Kelvin. The heating processes that generate, and sustain 
the hot corona have so far defied quantitative understanding, despite efforts spanning more 
than half a century ( |Kuperus et al.ligSlllGomezllQQOtlzrrkerliggSlpfmanpOOSIIKlimchukl 
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2006] |Taroyan & Erdelyi(2009) l. Efforts to establish the causes of solar atmospheric heating 
have produced a number of theories, of which two classes are most promising. The first 
holds for current dissipation, following reconnection events occurring throughout the at- 
mosphere in the form of micro-flare and nano-flare activity (Parker 1988 ; Priest & Schrijver] 
|1999 ; Fujimoto et al .|2011[ l. The second class predicts that the heating is driven by magneto- 
hydrodynamic (MHD) waves ( |Alfven| 1947HOsterbrockfl96Tl|Ionson 1 1 978[|HoIIweg| 199 1| >, 
most likely propagating from the lower atmosphere l |Roberts|2000 1. Turbulence can transfer 
energy from large scales to small scales. There is growing evidence that MHD turbulence 
and small-scale reconnections are explicitly linked, with current sheets acting as the dissipa- 
tive structures. Magnetic reconnection is also observed in simulations of MHD turbulence 
(Matthaeus & Velli 2011, Rappazzo et al. 2008). Two and three dimensional simulations 
show that magnetic structures tend to self-organize in space and time, creating a number 
of current sheets at their boundaries. The intermittency of microflares and nanoflares, com- 
bined with the small spatial scales over which the energy release occurs, allows us to refer 
to them collectively as turbulence (Einau di et al.| 199 6 ; Georg oulis et al.|1998| >. 

Kinetic processes and models that are based on wave particle interactions can not explain 
yet solar phenomena in a global scale, but are clearly applicable to small spatial scales. 
Kinetic processes and instabilities are highly relevant in models of plasma heating and solar 
wind acceleration. Hybrid models consider the acceleration of heavy ions by a spectrum 
of waves in the solar wind as well the influence of the ions on the wave structure and can 
reproduce the anisotropic velocity distributions of heavy ions ( |Scudder| 1 99"3a]b] |Matthaeusl 
|& Velli|201 lHOfman & Vinas|2007| i. 

A strongly related, and also unsolved question, involves the generation and acceleration 
of the solar wind, which often violently interacts with the magnetosphere and Earth's upper 
atmosphere. Again, there are strong arguments that MHD waves play a key role in contribut- 
ing to these important questions of space physics ( [Belcher| 1 97 1 1 1 Alazraki & Couturierj 1971 [ 
Gekelman 1999~} |Tsurutani et al.|19 99). The compelling theory of Alfven waves, which is 
the basis and starting point of most of the above seminal work, was developed by Hannes 
Alfven (Alfven 1942), for which he received a Nobel Prize in 1970. 

Fast and slow mode magneto-acoustic waves cause compression and rarefaction of the 
plasma as they propagate, with the magnetic and gas pressures acting as the main restor- 
ing forces, although kink modes are highly incompressible and the dominating restoring 
forces are magnetic tension-driven. Alfven waves are incompressible, with magnetic ten- 
sion the only restoring force for linear perturbations. If the perturbations are non-linear in 
nature, magnetic tension is the dominant restoring force. Alfven waves will cause primarily 
Doppler shifts, or specific magnetic field perturbations detectable with spectro-polarimetry 
in selected magnetically sensitive lines. Magneto-acoustic waves can cause both intensity 
variations and Doppler shifts. If the oscillating structure is spatially unresolved, the above 
effects will result in a broader line profile. 

A study of the dynamic properties of the interplanetary medium, carried out by the 
Mariner 5 spacecraft in 1967, has shown large -amplitude Alfven waves in micro-scale struc- 
tures ( Colem an| 1 967[ [Unt i & Neugebau er|1968| >. The largest amplitudes were identified in 
high velocity streams, but were also detected in low velocity regions. It is suggested that 
these are the remnants of waves generated near the solar surface, and have managed to prop- 
agate into the heliosphere without significant damping. The non-detection of appreciable 
magneto-acoustic wave power would suggest that these waves may have dissipated their en- 
ergy in the chromosphere and corona (Belc heT& Davis|1971| l. The Mariner observations, 
combined with the polar observations carried out by Ulysses, allowed a comprehensive view 
of turbulence in the heliosphere. Measurements of the solar wind parameters in the ecliptic, 
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have shown that the Alfvenic turbulence observed in high velocity streams, evolves towards 
a Kolmogorov-type spectrum with velocity shear helping to drive the process. The polar 
observations reveal a similar evolution but at a lower pace ( [Bruno &~Ca rbone 2005). These 
observations have prompted the development of wave-driven models to explain wind accel- 
eration and heating (Velli 2003) (see |Ofm an (2010) for a recent review). 

The motivation for the search for Alfven waves in the corona can be traced back to 
|Alfven| |1947| . Observations of Alfven waves in the magnetosphere (e.g. Keiling 2009) 
have prompted the search for the source of these waves in the lower solar atmosphere. With 
the advent of state-of-the-art instruments on-board satellite missions and ground-based so- 
lar telescopes, Alfven and transversal kink wave modes have started to be studied in the 
solar photosphere, chromosphere and corona at an unprecedented level of detail jAschwan-| 



|den et al.|199"9"l|Nakariakov et al.|1999[|Ofman & Wang 2008; Mcint osh etal.|2011) . The 
primary objective of this paper is to give an unbiased, but unavoidably personal account of 
the recent observational reports of Alfven waves in the solar atmosphere. This work will 
be underpinned by the linear magneto-hydrodynamic (MHD) theory that describes these 
waves in (non-)uniform magnetic structures (e.g. slabs and cylinders), which are considered 
to be simple, yet accurate descriptions of the actual magnetic building blocks of the solar 
atmosphere. We point out that the accuracy of slabs, cylinders, and pure linear modes in 
the description of realistic coronal active region loops is often doubted. Improved 3D MHD 
computational studies of coronal loops, in more realistic bi-polar active region geometries, 
has been recently performed by McLaughlin & Ofma n] l|2008|>, [O fman (2009a), |Selwa &| 
|Ofman| ( [20T0l ), |Schmidt & Ofman| ( |2011| ), and |Selwa et al.|p011a|b| l. 



2 The Transition Region and Corona 

2. 1 Non-thermal Line Broadening 

The global torsional oscillation of a coronal loop, with an axis that is not parallel to the line 
of sight, can produce periodic variations in the observed spectral line widths (McClements 
|et al.|1991[|Erdelyi et al.|1998[|Zaqarashvili|2003| l. The first observational efforts to detect 
Alfven waves in transition region (TR) and coronal structures dates back to the 1970s, with 
observations coming from the SKYLAB and SMM missions. Analysis of UV spectra, ob- 
tained with the Naval Research Laboratory spectrograph, revealed non-thermal velocities in 
the ranee of 10-30 km s" 1 ( 
from these studies is that, if the excess line broadening is to be attributed to Alfven waves, 
these waves cannot account for the heating. McCleme nts et al.] I |1991| l suggested that un- 
ambiguous evidence for Alfven waves in the TR can be provided by tracking coronal loops 
as they move from disk center to the limb. The non-thermal widths of emission lines from 
these loops would be expected to increase as a function of distance from disk centre. These 
line widths should become maximum at the limb, where the magnetic field lines are pre- 
dominantly perpendicular to the observer's line of sight. Indeed, SUMER observations have 
shown a centre to limb increase in the non-thermal line widths of upper chromospheric and 
TR lines, while the excess broadening of coronal lines is only marginal (see e.g. Erdelyi et 
|al.| 19 98 ; Bane rjee et al.|2 007 and references therein). However, chromospheric and transi- 
tion region lines are known to suffer from opacity effects. Thus, the center-to-limb variation 
of the line widths can also be attributed to increased opacity broadening towards the limb. 

The line broadening could also be due to the effects of bulk outflow motions that cause 
differential Doppler shifts along the observer's line-of-sight. The height profile of the tem- 



Doschek et al. |l976| [Cheng et al. 1979). The main conclusion 
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perature distribution, bulk outflow motions, and the geometry of superradial expansion will 
enhance the Doppler broadened line profile. These effects are exacerbated by an anisotropic 
temperature distribution, particularly when T_l >> Tn 



|& SolanE] l |2006| > also find that at distances greater 



l |Akinari|2007[|Ofman|20T0a| l. 
that IKf?) from the solar suria 



Raoua i 

^0 ironi me soiar surface, the 
widths of some transition region and coronal lines depend on the adopted electron density 
stratification. This became apparent in polar coronal hole observations of O VI and Mg X 
lines. This may indicate that the need for anisotropic temperature/velocity distributions may 
not be as important as previously thought. 

In a follow up investigation, |Chae et al.| ( fT998| > argue that the small center-to-limb vari- 
ation can be due to network cell structures where the non-thermal velocities may be consid- 
ered isotropic. The authors concluded that if the small increase in the non-thermal velocities 
of EUV lines is due to Alfven waves, the limited temporal resolution of the observations 
would imply that the wave periods will be shorter than a few tens of seconds. They used 
the non-thermal velocities to derive a wave flux, which was then compared with the radia- 
tive losses in the 10 4 - 2xl0 6 K temperature range. Their comparison reveals that there 
is sufficient Alfven wave flux to balance the radiative losses for temperatures in excess of 
3x 10 4 K. For temperatures lower than 3x 10 4 K, the wave energy is significantly smaller 
than the radiative losses. However, this is inconsistent with Alfven wave flux which is ex- 
pected to remain constant (or decrease) as the temperature increases ( |Chae et al.|1 998 l. This 
rapid increase in energy as a function of temperature is a direct consequence of the small 
increase in the non-thermal velocities found in this temperature range (Fig. [TJ. 
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2.2 The CoMP View of the Corona 



Direc t measurements o f the coronal magnetic field have been uncertain, and notoriously dif- 



ficult. Lin et al. 1 2000 1 have used the near infrared magnetically sensitive Fe XIII 10747 A 
line to measure the Stokes V profiles above active regions, where they determined longitu- 
dinal magnetic field strengths in the range of 10 - 33 G. 

The state-of-the-art Coronal Multi-channel Polarimeter (CoMP; Darnell et al. 2003} has 
been employed behind a coronagraph at the Hill-top facility at the National So lar Observa 



tory in New Mexico, to provide the states of polarization in the 10747 A line I Tomczyk et 
|al.|2007| >. The observing setup allows measurements of spectral line intensities, line-of-sight 
velocities, and line widths in the corona between 1.05 and 1.35 solar radii. The observations 
show clear evidence for line-of-sight velocity variations with a period of approximately 5 
min (3.5 mHz). Simultaneous measurements of linear polarization have revealed the plane- 
of-the-sky azimuthal component of the coronal magnetic field. jTomczyk et aT7| ( |2007} l sug- 
gested that these oscillations are the signatures of Alfven waves in the corona for the follow- 
ing reasons: (a) The waves have phase speeds of about 2000 kms -1 , which is considerably 
higher than the coronal sound speed, (b) The waves propagate along density structures that 
are thought to follow the magnetic field lines, and (c) There is no evidence for appreciable 
intensity oscillations which indicates that the detected wave modes are incompressible. The 
lack of intensity oscillations may be due to the low spatial resolution that would average out 
any small scale variations (see later). 

The Fourier spectrum of the Doppler velocities associated with these waves points to 
a photospheric driver (i.e. p-modes). However, since their frequency is below the chromo- 
spheric cut-off, the only way that they can propagate upwards into the corona is through 
inclined magnetic structures l |De Pontieu et a l. 2004; Erdelyi 2006a b|. They would also 
need to be converted from pressure-driven modes into non-compressible transversal Alfven 
modes. 

The interpretation of the CoMP results, in terms of Alfven waves, has been criticised by 
e.g. | Van Doorsselaere et al.| ( |2008| ), who argue that these observations show the collective 
motion of the plasma from its equilibrium position, thus causing a periodic variation in 
the line-of-sight velocity. The authors suggest the best interpretation is the fast kink mode. 
We emphasise that even a small intensity perturbation would also point towards the fast 
kink MHD mode. However, one has to consider how the scales of the photospheric driver 
compare with the scales of individual threads. If the photospheric driver has dimensions 
smaller than the threads, the idea of kink waves is favourable. As we are dealing with waves 
in chromospheric and coronal structures, the kink interpretation becomes more meaningful 
if the structures are stable over several wave crossing times. 

Irrespective of the actual interpretation, there is no doubt that the observational findings 
of Tomczyk et al. (2007 ) are sound, and clearly demonstrate the ubiquitous nature of waves 
in the solar corona. In a follow up investigation, Tomczyk & Mcintosh ( 2009 ) performed 
time-distance seismology of the solar corona, and constructed a k-co diagram of the region. 
The power spectrum of the coronal velocity perturbations shows an excellent agreement 
with that of the solar p-modes, which supports the suggestions that these waves can be 
transported into the corona along inclined flux tubes originating in the photosphere (Fig. [2}. 
However, the high-frequency spectrum is distinctly different from the photospheric power 
spectrum, and indicates that the high-frequency component of these waves may be generated 
locally. An estimate of the wave energy reveals that it is significantly smaller than the energy 
required to balance the coronal radiative losses. However, it should be emphasised that the 
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Fig. 2 The Fourier spectrum of the Doppler velocities in the corona (black line) and the photospheric p-mode 
spectrum (green line). Power-law fits to the low (red) and high frequency (blue) parts of the coronal spectrum 
are also shown (reproduced from Tomczyk & Mcintosh 2009). 



low spatial resolution of CoMP, limited to 4.5 arcseconds/pixel, overlooks small-scale solar 
structures, and therefore significantly underestimates the energy carried by these waves. 

The higher spatial resolution of the Atmospheric Imaging Assembly on board the Solar 
Dynamics Observatory (SDO; Boerner e t al.|20lT| has indeed shown transverse oscillations 
in the corona with amplitudes of around 20 kms -1 . The amplitude of these motions is a 
factor of 40 higher that those found by CoMP, and may be sufficiently energetic to balance 
the radiative losses from the quiet solar corona ( |McIntosh et al.|2011| >. 



2.3 Flare induced Alfven Waves 

Solar flares often induce temperature rises in the chromosphere of the order of 100 - 200 K. 
This requires a significant energy input, amounting to 10 - 20 ergcm~ 3 s^ 1 at these at- 
mospheric heights. The dense nature of the lower solar atmosphere, when compared to the 
corona, makes it difficult to reconcile this amount of heating from an individual energy 
release in the corona (Mach allo et al.| [T978). This is often called the "number problem", 
whereby the high total number of electrons required to be accelerated down into the lower 
atmosphere do not exist at coronal heights ( |Fletcher & Hudson |2008| >. 

An alternative interpretation to explain lower-atmospheric heating is that of flare-induced 
dissipation of Alfven waves. Emslie & Sturrock ( 1982) have suggested how Alfven-wave 
perturbations may be generated along current sheets, linking the corona to the lower solar 
atmosphere, in the aftermath of an impulsive flare event. As the group velocity of Alfven 
waves will be parallel to the magnetic field, the energy contained within the mode will re- 
main trapped, providing the wavelength is small compared with the scale length of field 
and plasma variations. Thus, energy can be readily transported down into the lower solar 
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Fig. 3 The coronal magnetic field launches a torsional Alfven wave pulse into the lower atmosphere, as well 
as a fast-mode wave pulse. That fraction of the Alfven wave energy that survives into the chromosphere can 
also lead to stochastic acceleration there. The wave will be partially reflected from the steep gradients in the 
chromosphere (not shown) and re-enter the corona (reproduced from Fletcher & Hudson 2008 1. 



atmosphere, where dissipation mechanisms convert the energy into localized heating. The 
Poynting flux of the waves is proportional to the Alfven speed and the square of the magnetic 
field perturbation. X-class flare energy requirements would suggest a Poynting flux of ap- 



proximately 10 11 ergcm~ 2 s _1 . The magnetic field variations, as high as 200 Gauss I Sudol 
|& Harvey|2005| l, observed during the most energetic solar flares, combined with the Alfven 
speed in the corona, implies that there is sufficient energy to meet the flare requirements 
(Fletcher & Hudson|2008l l. 

What triggers Alfven waves during flare events is still an open question. However, cur- 
rent understanding suggests they can be generated during the reconnection process, when the 
magnetic field relaxes from its stressed pre-flare state (Aulanier et al.|200 6). In this regime, 
magnetic field lines may still be highly distorted from a potential configuration, providing 
the necessary high-tension force required to induce Alfven waves (Fig[3j Fle tcher & Hudson] 
|2008) . 

A high degree of Alfven wave interaction will occur as these waves impact the lower 
solar atmosphere. Goedbloed & Halberstadt ( 1994) have shown how incident Alfven waves 
will produce a fast- wave spectrum containing rapidly varying amplitudes. It is these waves 
which can be locally damped through ion-neutral or resistive damping, or form a turbulent 
cascade which propagates deeper into the photosphere (Emslie & Mach ado|1979| l. 

Magnetic reconnection can also produce small scale X-ray jets when a burst of solar 
plasma is driven into the corona at high velocities (Shibata e t al. 1992). These X-ray jets, 
first observed by Yohkoh, have been studied in greater detail by the Hinode X-ray Telescope 
(XRT). The XRT studies ident ified a large numbe r of X-ray jets moving with velocities in 



the range of 200 - 800 km s 1 i Cirtain et al. 2007 1. The high end of these values is similar to 



the Alfven velocity in the corona. The suggestion that an Alfven wave may be generated in 
these reconnection sites, prompted |Cirtain et al.|(200"7| to search for transverse oscillations 
associated with the location of the jets. They identify transverse motions with a period of 
approximately 200 s, which seem to be related to the high velocity component of the outflow. 
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3 The Photosphere and Chromosphere 

3.1 Transverse Alfven Waves in Spicules and Fibrils 

The wealth of observations from a wide range of ground-based and space-borne missions 
have shown that the solar chromosphere is permeated by spicules and fibrils. Spicules are 
dynamic, straw-like magnetic structures that are best identified in the Ca II H & K lines, 



Na D3, Ha and Ca II 8542 A, and can be divided into two distinct classes ( Zaqarashvili & 
Erdelyi 2009). Many detections of transversal waves in spicules have been tied to Type I 
spicules, which are ubiquitous throughout the solar atmosphere, and longer lived (^10 min) 
than their jet-like Type II counterparts ( |De Po ntieu et al. 2007). However, the classification 
of spicules into two separate types has been recently called into doubt (Zh ang et al.|2012) . 
Formed in low plasma f3 regions, the axis of spicules and fibrils outlines the direction of 
the magnetic field. Transversal oscillation in these structures, interpreted as Alfven waves, 



|2008| l. These authors used the Solar Optical Telescope (SOT; [Ts"uneta et a! 



corona loops 


Ofman & Wang 


, Tsuneta et al. 


2008)onHinode 



to analyze a time series of Ca II H (3968 A) images. The Hinode observations clearly show 
that spicules undergo transverse motions, with periods in the range of 100 - 500 s and 
amplitudes of 10 - 25 kms -1 . The oscillations are evident in traditional x — t plots, where 
the intensity is plotted as a function of time in a direction perpendicular to the axis of the 
spicule. The estimated periods can be uncertain as they are comparable to the lifetime of 
the spicules. Estimates of the energy flux carried by these waves reveals that it is sufficient 
to balance the radiative losses from the corona and power the solar wind. Recently, |He et| 
|aT| ([2009 ) used a high-cadence mode on Hinode to identify high-frequency oscillations in 
spicules with periodicities less than 50 s. These waves are interpreted in terms of upwardly 
propagating Alfven waves with phase speeds >75 km s , and wavelengths less than 8 Mm 
l |Martinez et al.|2009) . 

While transverse kink motion has been shown to be abundant in spicule observations, 
the underlying cause of the periodic motions has remained speculative. Overshooting of con- 
vective motions in the photosphere, granular buffeting, rebound shocks, and global p-mode 
oscillations have all been suggested as candidates for the creation of spicule oscillations 
( |Roberts|[T979l |Sterling & Hollweg|[T988l |Vranjes et aT][2008l i. |De Pontieu eTaT] pOO^ 
have suggested how photospheric p-modes can develop into shocks, and drive the dynamics 
of chromospheric spicules. The authors indicate how, even though low frequency p-modes 
become evanescent in the upper photosphere, these waves can propagate upwards into the 
upper chromosphere and corona through inclined flux tubes. Flux tubes with high inclina- 
tion angles enhance the leakage of p-modes into the upper layers of the solar atmosphere. 
Recently, |Jess e t al. (2012) have presented unambiguous evidence of how the generation of 
transverse oscillations in spicule structures is a direct result of the mode conversion of longi- 
tudinal waves in the lower solar atmosphere. Furthermore, these authors show how spicules 
are intrinsically linked to photospheric flux concentrations, with their spectrum undoubtedly 
including signals from both the 3-minute and 5-minute p-mode oscillations (Fig.[4]i. 

The question of how these waves dissipate remains open. Wave energy damping in the 



Hollweg & Yang|1988||Poedts et al. 
1994||1995||Ofman & Davila|1995 ; 


1989||1990[|Steinorfson & Davila|1993||Ofman et al. 
Erdelyi and Goossens 1995; Ruderman and Roberts 


2002 ; Goossens et al.|2011[l, phase mixing 


(Heyvaerts & Priest 1983; Ofman & Aschwan- 


den 2002 ), Alfven turbulence ( Cranmer & van Ballegooigen 2005 ; van Ballegooijen et al. 
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Fig. 4 Simultaneous images of the Ha core (chromosphere; upper left) and G-band (photosphere; upper 
right), where dashed white lines in the Ha core image highlight the spatial positions where x — t cuts were 
made. The blue trace in the G-band image denotes the position of a chromospheric spicule (visible in the 
left panel as a dark, straw-like structure), where one end is anchored into the photosphere above a MBP. A 
red contour indicates the location where a high concentration of longitudinal oscillatory power is present. 
The lower panel is a sample Ha time-distance cut, obtained 4000 km (5.5 arcsec) from the underlying MBP, 
revealing an abundance of periodic transverse motions in the solar chromosphere. White arrows highlight a 
trough and a peak of a typical transverse oscillation. The scale is in heliocentric coordinates, where 1 arcsec 
approximately equals 725 km (reproduced from Jess et al. 2 012) . 



201 1 ) and/or mode conversion ^Suzuki & Inutsuka|2005| l. |Velli et al.| ( fT989] > suggested wave- 
reflection in the steep density gradients of interplanetary medium as a source of the turbulent 
cascade of unidirectional Alfven waves. This idea was applied to the lower solar atmosphere 
by Mattha eus et al.| fl999 ) who proposed that the Alfven waves generated in the photosphere 
and chromosphere interact with their reflections to drive the effects of turbulence offering a 
viable mechanism for heating the corona in the open field regions. Motivated by the CoMP 
observations, Verdini et al. (2009 ) have re-examined the turbulent spectrum from the base of 
the corona and out to 17 R©. They find that turbulent dissipation can account for at least half 
of the heating required to sustain the solar wind. Reflection driven turbulence can therefore 
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play a vital role in the acceleration of the fast solar wind. However, the main limitation of 
these models is the use of incompressible reduced MHD equations. Thus, the likely impor- 
tant coupling between Alfven waves and compressive modes is neglected in these models. 
|He et al.| (2009) suggest that Alfven waves could instead be damped by cyclotron resonance 
in the heliosphere, out to 60 R©, a mechanism that is also favourable for the heating of the 
solar wind (see reviews by Marsch ( 2006} ; [Ofman| p0T0"| l) . 

The evolution and dissipation of Alfvenic perturbations in force-free magnetic struc- 
tures has been studied in detail by Malara et al. (2000), who examined the competing effects 
of phase mixing and three dimensional effects (i.e. exponential separation of field lines) in 
wave packet propagation. They find that in a chaotic magnetic field configuration, the dis- 
sipation is dominated by three dimensional effects, while in simpler equilibrium situations 
both phenomena contribute in spatially separated regions. We point out that this approach 
is limited by the assumption of a incompressibility. Large amplitude Alfven waves are sub- 
ject to parametric decay as a result of their coupling with acoustic and transverse waves. 
The energy of an Alfven wave propagating in a turbulent medium can be transferred to 
other wave modes. Malara & Velli ( 1996) have shown that the parametric instability of non- 
monochromomatic waves decreases quadratically with departures from the monochromatic 
case. Numerical simulations of the time evolution of large amplitude Alfven waves have 
confirmed the development of this instability in the solar wind (Malara et al. 2000a; Del 
|Zanna et al.|200T] ). 

Spicule oscillations may also be attributed to the transversal kink MHD mode ( |Za-| 
qarashvili & Erdelyi 2009 ). The observational signatures of the guided-kink magneto-acoustic 
mode are the swaying, transversal motions of magnetic flux tubes, similar to the oscillatory 
behaviour observed in spicules (Erdely i & Fedun|2007| ). [De Pontieu et al.| ( [2007[ ) have ar- 
gued that, due to their short lifetime, the spicules cannot act as stable wave guides and the 
kink mode cannot satisfactory explain the Hinode observations. This apparent contradiction 
may be resolved by considering the spicule as sliding along the oscillating field lines. 

In view of the recent Hinode SOT analysis of spicules, Mci ntosh et al.| ( |2~008[ ) have re- 
examined the SUMER results of transition region line widths by combining the results of 
forward modelling with spicule observations. The authors conclude that the C IV intensities 
and non-thermal velocities can be explained by the superposition of longitudinal mass flows 
and 'Alfvenic motions". It remains unclear what is meant by 'Alfvenic motions", since the 
eigensolutions of MHD equations are not labelled commonly as such in theoretical studies. 
On the other hand it is not evident that linear MHD eigenmodes of cylindrical magnetic 
structures can well-describe all realistic coronal conditions. The apparent isotropy in the 
center-to-limb variations may be explained by the similar magnitude of observed longitudi- 
nal and transverse velocity components. 

Solar prominences, the heavy and cool elongated magnetic structures supported by 
mainly horizontal magnetic fields in the solar corona, can support a wide variety of MHD 
waves. Analysing the properties of these waves can give an unprecedented insight into the 
small-scale structuring of these features, with clues about the magnetic field strength through 
means of solar magneto-seismology (Nakariak ov & Verwichte|2005| >. 

Prominence oscillations cover a wide range of amplitudes. The large amplitude velocity 
oscillations, >20 kmr 1 , are thought to be excited by Moreton waves produced by flares 
(Moreton & Ramsey 1960). These oscillations are observed in both the transverse and lon- 
gitudinal directions, with periods in the range of tens of minutes to a few hours (Trip athi etj 
|al.|20 09l. The potential of large-amplitude prominence oscillations for deriving the vertical 
magnetic field was first explored by Hyder ( 1966), who obtained field strengths of tens of 
Gauss. More recently, Vrsna k~et al.| ( [2007] ) used prominence seismology techniques to in- 
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fer the Alfven speed, as well as the azimuthal and axial components of the magnetic field 
strength. The small amplitude oscillations, <3 kms -1 , are usually restricted to only part 
of the prominence, and are thought to be excited by the photospheric and chromospheric 
oscillations. 

The vertical oscillatory motions observed in prominences have been interpreted as trans- 
verse waves occurring on the horizontal magnetic field lines of these structures. The oscil- 
lations have periods of 2 - 4 min, and are in phase along the entire length of the prominence 
thread, as observed in the plane of the sky (O kamoto et al.|2007] >. The lack of any additional 
information make the exact wave mode difficult to disentangle. We therefore believe that 
the suggestion that these oscillations are due to Alfven waves may be somewhat premature. 
This is further supported by numerical simulations by e.g. Erdelyi & Fedun (2007). For a 



recent review on prominence oscillations see Arr egui et al . 1 1 20 1 2 ) 



3.2 Torsional Alfven Waves 

The periodic motions of photospheric footpoints, in an axially symmetric system, can excite 
torsional (m = 0) Alfven waves (see 5j5]of this review and Ruderman et al. 1997; Ruderman 
1999). The observational signature of torsional Alfven waves, with a velocity component 
along the observers line-of-sight, will arise from torsional velocities on small spatial scales. 
It is suggested that these torsional velocities will produce an observable full-width half- 
maximum (FWHM) oscillation in a given line profile. At the linear limit of infinitesimally 
small amplitude, Alfven waves are incompressible, and should therefore exhibit no periodic 
intensity variations. 

The detection of such torsional motions in the lower solar atmosphere has, until recently, 
remained elusive. |Jess et ak| |2009 ) used the Swedish Solar Telescope, equipped with the 
Solar Optical Universal Polarimeter (SOUP; Title 1984), to carry out imaging spectroscopy 
of the lower solar atmosphere. The tuneable nature of SOUP allowed complete sampling 
of the Ha profile, and multi-object multi-frame blind deconvolution ( |van Noort et al.|2 005 1 
image reconstruction was implemented to remove small-scale atmospheric fluctuations. 

The authors investigated a large conglomeration of photospheric magnetic bright points 
(MBPs), covering an area of 430,000 km 2 . FWHM oscillations of the Ha line profile, in 
addition to line-of-sight Doppler blue-shift velocities of 23 kms -1 , provided evidence for 
the presence of torsional (m=0) Alfven waves. The progression of the Doppler-compensated 
Ha FWHM is shown as a function of time in Figure [5] with a schematic diagram of the 
observed wave mode shown in Figure [6] A lack of periodic intensity fluctuations in their 
dataset further strengthened the authors conclusions. 

Cauzzi et al.|(|2009l have used the Interferometric Bi-dimensional Spectrometer (IBIS; 



|Cavallini|2006 1 to perform imaging spectroscopy of the solar chromospheric network in the 



Ha and Ca II 8542 A line. They find that an increased Ha line width correlates well with 
the Ca II minimum intensity. This is interpreted as evidence for shock heating in the mag- 
netic network which raises the possibility that acoustic shocks can introduce an oscillatory 
behaviour in the observed Ha line width. To investigate this possibility further, |Jess et ak] 
( 2009 ) carried out a number of additional tests to help verify or refute their original conclu- 
sions. Since torsional Alfven waves will produce a FWHM oscillation that is 180° out of 
phase at opposite boundaries of the waveguide ( |Van Do orsselaere et al. 2008 ), the authors 
investigated the instantaneous Ha profile at opposite edges of the MBP group. Through 
comparison to an at-rest profile, the nature of wing broadening is clearly visible in the lower 
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Fig. 5 A wavelength versus time plot of the Ha profile showing the variation of FWHM line width as a 
function of time. The arrows indicate the positions of maximum amplitude of a 420-s periodicity associated 
with the MBP group (reproduced from Jess et al. 2009 1. 



panel of Figure [7] These findings reinforce the conclusion of | Jess et al.| ( |2009] > that the ob- 
served FWHM oscillations are due to the presence of torsional Alfven waves. 

The energy flux carried by the Alfven waves can be estimated using the WKB approxi- 
mation, 

Fa = Pv 2 va , (1) 
where p is the mass density of the flux- tube, v is the observed velocity amplitude, a nd va the 
Alfven speed. For a mass density of p w 1 x 10 -6 kg m -3 ( Vernazza et al. 1981 1, and uti- 
lizing the derived wave parameters described by |Jess et al.| l |2009| >, the localised energy flux 
in the chromosphere is ^150,000 Wm~ 2 . It is estimated that, at any given time, wl.6% of 
the solar surface is covered by bright point groups (Sanchez Almeida et al. 2010), producing 
a global average of 2400 Wm~ 2 . This value is significantly above the energy required to 
heat the corona and power the solar wind. We emphasise that the above energy estimate is 
correct in the absence of vertical gradients or when they can be neglected in the WKB limit. 
The presence of a soalr wind enhances the transmission of low frequency waves, while high 
frequency waves have a transmission similar to static models (Heine mann & Olbert[1 980; 
|Velli|1993[ ). 



4 Wave Drivers 



Turbulent convection leads to the stochastic excitation of p-mode oscillations, which are 
ubiquitous in the lower solar photosphere, and can be detected in both spectral line intensi- 
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Fig. 6 Schematic diagram of an expanding magnetic flux tube, sandwiched between photospheric and chro- 
mospheric intensity images, undergoing a torsional Alfven perturbation that propagates in the vertical direc- 
tion. At a given position along the flux tube, the Alfven displacements are torsional oscillations that remain 
perpendicular to the direction of propagation (not to scale; reproduced from Jess et al. 2009 1. A divergent 
magnetic geometry leads to wave amplitude increase with height and coupling to compressional waves (Of-| 
|man&Davila|1998> . 

ties and Doppler shifts. The frequency of p-mode oscillations are in excess of 1 mHz, with 
the strongest power detected at ^3 mHz (5 min). However, the p-modes can be modified 
by the magnetic field (for reviews with references on this subject see e.g. Erdelyi 2006a; 
Pinter & Erdelyi 201 1 ). In magnetic regions such as sunspots and pores, the amplitudes of 
the 5-min oscillations are suppressed, while high frequency oscillations, above the acoustic 
cut-off, have enhanced power. The suppression of low frequency oscillation occurs on very 
small scales (a few arcseconds), while the enhanced high frequency power occurs primarily 
in patchy halo-like areas surrounding the strongest magnetic fields (Thomas & Stanchfield 
[20D0l >. 

The solar magnetic network, defined by the boundaries of super-granular flows, shows 
evidence for MHD waves with frequencies in the range of 1 - 4 mHz. The internetwork 
has frequencies in excess of 4 mHz, which have been modelled as acoustic waves driven 
by photospheric velocities (Carlsson & Stein 1997). Network bright points are formed at 
the boundaries of super-granules, and show a strong correspondence with the underlying 
magnetic field. The granular buffeting of photospheric flux tubes can create transverse waves 
which propagate along the field lines and into the upper atmosphere. The models of Musielak 
& Ulmschneider ( 2003 1 suggest that the low frequency oscillations in the magnetic network 
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Fig. 7 Top Panel: Simultaneous plots of the Ha line profile at opposite edges of the MBP group. The red line 
shows a profile which is heavily broadened in the red wing of the Ha profile, while the blue line displays a 
large degree of non-thermal broadening in the blue wing of the Ho profile, at the opposite edge of the MBP 
group. A rest Ho profile is displayed as a solid black line, with the Ha core rest wavelength plotted as a 
vertical dashed black line. Lower Panel: The same broadened profiles plotted above, but now displayed as 
intensities relative to the rest Ha profile. It can be seen that the resulting profiles are anti- symmetric about 
the Ha core rest wavelength of 6562. 808A, as indicated by the vertical dashed black line. This out-of-phase 
phenomena, at opposite edges of a flux tube undergoing torsional displacements, is a characteristic of the 
m=0 Alfven mode. 



are due to transverse waves that are generated by the granular buffeting of photospheric flux 
tubes. Density stratification will lead to an increase in the wave speed as they travel into the 
upper photosphere and chromosphere, where mode transformation can occur. Transverse 
waves can transfer power to longitudinal modes, which then develop into shocks and heat 
the surrounding plasma. A number of criteria may be used to identify if mode transformation 
has occurred within a given network bright point (Kalkofen 1997; Hasan & Kalkofen 1999; 
Vigee shet al.|201 1) . These include: 

- Power at a particular frequency, v, as identified in the upper-photosphere/lower-chromosphere, 
should also be identified in the upper chromosphere, but with reduced amplitude. 

- Power at double frequency, 2u, should also appear in the upper chromosphere. 
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- Waves remaining at, or above the transverse cut-off, should propagate at approximately 
the sound speed. 

- The oscillatory signal should be quasi-periodic, as a result of the granular buffeting. 

The studies of |McAteer et al.| < |2003] > and |Bloomfield et al.|(2004) have provided the ob- 
servational evidence for mode coupling in such magnetic networks. In these examples, the 
authors provide evidence for transverse waves with frequencies in the range 1 - 2.5 mHz, 
transferring power to longitudinal waves with frequencies of 2 - 5 mHz. We emphasise that 
the above observational signatures could also appear as a result of nonlinearity. 

It has been demonstrated, both analytically and theoretically, that when the magnetic 
pressure is approximately equal to the gas pressure (plasma /3 = 1) in a solar flux tube, 
longitudinal to transverse mode coupling may also occur {be Moortel et al. 2004). Until 
recently, this form of "reversed mode coupling" had not been verified observationally. How- 
ever, [jesFetan]([20T2| have suggested that longitudinal- to-transverse mode conversion may 
account for the transversal waves which are ubiquitous in the solar chromosphere. Jess et al.] 
(2012) utilised high spatial and temporal resolution observations, obtained with the Rapid 
Oscillations in the Solar Atmosphere (ROSA; |Jess et aL||2010| > imager at the Dunn Solar 
Telescope, to investigate the origin of transverse oscillations in chromospheric spicules. The 
authors present strong evidence that transverse oscillations in Type I spicules, with periods 
in the range of 65 - 220 s, originate in longitudinal oscillations of photospheric magnetic 
bright points, at twice the transverse period. It must be emphasised that this form of mode 
coupling occurs in the reverse direction of that reported by |McAteer et al.| < |2003[ >, i.e. lon- 
gitudinal to transverse. The mechanism for the mode transformation revolves around the 
nature of the photospheric driver, whereby a 90° phase shift across the body of the spicule 
creates gradients in pressure, which displace the spicule axis, generating the observed kink 
modes. These periodic pressure differences also induce a frequency doubling of the coupled 
wave mode, in addition to producing periodic compressions and rarefactions along the body 
of the spicule (i.e. the sausage mode). 

It is known that convective motions, intrinsic to photospheric granulation, transport an- 
gular momentum and agitate plasma in the outer layers of the Sun's atmosphere. Solar ma- 
terial rises in the granules and sinks in the inter-granular lanes. The plasma that returns to 
the solar interior has increased angular momentum, and therefore generates an increased 
amount of vorticity at the edges of the lanes (Stein & Nordlund 1998). The idea of vortices 



as the source of torsional wave modes has been proposed by Hollweg et al. ( 1982) and |Velli| 
|& Lie wer (1999). Small scale vortex motions, with dimensions less than «500 km, have 
recently been observed in photospheric G-band images of MBPs (Fig.[8j B onet et al.|2 008 1. 
These small scale vortices have lifetimes of a few minutes, and appear to trace the pattern 
of a super-granular cell. They seem to be formed when two MBPs move towards a common 
point, and rotate relative to one another. Wedemeyer-Bohm & Rouppe van der Voort ( 2009 1 
have observed similar vortices in the chromospheric Ca II 8542 A line core. The location of 
the chromospheric vortices coincides with photospheric MBPs, although there is no current 
evidence for simultaneous swirl events in the corresponding photospheric locations. 

Recent radiative MHD simulations have shown photospheric vorticities in intergranular 
magnetic field concentrations ( Shelyag et al.|2011[ ). The connection between MBP vortices 
and upper atmospheric swirls appears to be magnetic in nature. Simulated G-band images 
show a direct comparison between magnetic vortices and rotary motions of photospheric 
bright points. These swirly motions of plasma in the intergranular magnetic field concentra- 
tions could be responsible for the generation of different types of MHD wave modes such 
as kink, sausage and torsional Alfven waves. 
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Fig. 8 A sequence of G-band snapshots showing a photospheric vortex. Two pairs of MBPs are marked as 
squares and crosses to help the reader identify the rotational motions. The snapshots have been taken with a 
cadence of 15 s. The bar in the lower left corner corresponds to 1,000 km on the solar surface (reproduced 
from |Bonet et al.|2008) . 



5 Theory 

5.1 Alfven Waves in Uniformly Magnetised Plasma 

Within the linear limit of infinitesimally small wave amplitudes in uniform, finite beta mag- 
netized plasma, there are three distinct types of MHD waves that can be mathematically 
described as: slow and fast magneto-acoustic waves, and Alfven waves. The first two types 
of wave have an acoustic character modified by the magnetic field, whereas the Alfven 
waves exist purely because of the presence of a magnetic field. The three physically distinct 
MHD wave modes exist due to the following reasons. The MHD description of a magnetised 
plasma is using macroscopic quantities, let's say density (p), velocity (v), magnetic field 
(B) and a thermodynamics variable, say pressure (p). In ideal MHD (i.e. no dissipation or 
non-adiabatic processes) these quantities are coupled through a set of eight non-linear par- 
tial differential equations (PDEs) governing mass conservation, momentum conservation, 
governing magnetic flux and the energy conservation. If this is a closed system, in a station- 
ary state one may expect eight eigenvalues corresponding to eight different physical wave 
modes. However, there is the solenoidal condition, that will reduce the number of indepen- 
dent PDEs to seven, resulting in seven eigenvalues. One eigenvalue, u = 0, turns out to be 
describing the entropy wave that does not carry information, and may be disregarded, leav- 
ing six possible independent eigenvalues corresponding to six further independent modes. 
An interesting intrinsic feature of the remaining eigenvalues is that they appear squared, i.e. 
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there is no distinction between forward and backward propagating MHD waves in a static 
(i.e. vo = 0) magnetised plasma. Hence, there are three physically distinct MHD waves. 
These waves are a complete set and their appropriate linear superposition can be combined 
into any linear MHD wave perturbations. This also means that the individual eigenmodes 
(slow, fast magneto-acoustic and Alfven waves) are orthogonal to each other and there is 
no linear transformation that can transform one MHD eigenmode into the other. This latter 
property is very important when identifying MHD modes in observations. Often, an accu- 
rate mode identification is possible only by considering very carefully the mode dispersion 
or the phase relation between the perturbations involved. An interesting and educational ex- 
ample of this approach is e.g. by Fujim ura & Tsune"ta| ( |2009| . The problem becomes even 
more subtle in structured plasmas (e.g. magnetic cylinder or slab), where the geometry of 
the problem introduces an infinite set of MHD eigenmodes. 

5.2 Alfven Waves in Structured, Non-uniform, Magnetised Plasma 

A natural way to improve the simplistic view of a uniform MHD plasma (c.f. j ]5.1| >, in order 
to narrow the gap between theory and the observed wave processes taking place in the solar 
atmosphere, is to add non-uniformity. One way to achieve this goal is to introduce struc- 
turing. This may be done in many ways, but both the slab and cylindrical geometry seem 
to be a popular approach to model the building blocks of the solar atmosphere. We empha- 
size that the conservation of distinct modes in the real coronal plasma is not guaranteed, 
since the coronal plasma is nearly always more complex in nature than a collection of slabs 
and cylinders (Ofman 20091. Another possibility to introduce non-uniformity is to consider 
inhomogeneity. Inhomogeneity may be either along the waveguide (for example for a verti- 
cal magnetic waveguide this could be the gravitational or magnetic stratification) or across 
(e.g. for a vertical flux tube, this may mean inhomogeneity in the radial direction). It was 
always expected that once MHD waves are generated, they will easily propagate along mag- 
netic flux tubes (homogeneous or non-uniform), the building block of the solar atmosphere 
(Fig.[9^), or along magnetic field lines at constant magnetic surfaces (Fig.[9j>). 

MHD wave theory has been around for a while, in particular for linear magneto-acoustic 
and Alfven waves in structured magnetised plasmas (cf., |Roberts[198 1, Edwin & Roberts] 
|1982| for slab and Edwin & Roberts 1983 for cylindrical geometries). The propagation of tor- 
sional Alfven waves in vertical magnetic tubes have been studied by, e.g., Hollweg ( 1981 1; 
Holl weg et aT] ( [T982]> a nd more recently in the context of filtering torsional Alfven waves 
by Fedun et al. (201 lb). In order to overview briefly our current understanding of what kind 
of linear MHD waves are supported by MHD structures, we will assume a simple cylindri- 
cal geometry with some non-uniform radial inhomogeneity. We would also point out that 
it would not be practical to single out the Alfven wave mode, as there is coupling between 
the linear MHD modes. The waves have to be investigated in a co-existing context. In na- 
ture, however, there may be a physical situation where the excitation mechanism (or driver) 
would favour one of the modes. Moreover, other wave modes will be produced due to linear 
(in non-uniform media) and non-linear (even in an initial inhomogeneous plasma) coupling. 
In order to have the problem tractable for the purpose of wave mode identification, with a 
focus on Alfven wave, we assume a uniform magnetic twist here and will follow the key 
steps of |Erdelyi & Fedun| ( |20T0^ . 

Let us consider a magnetically twisted straight flux tube embedded within a uniformly 
magnetised plasma environment in cylindrical geometry. We note, such structures are al- 
ways unstable, unless they are supported by appropriate boundaries and the twist does not 
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Fig. 9 (A) Magnetic flux tube showing two snapshots (at positions Z\ and 22 ) of Alfven wave perturbations 
propagating in the longitudinal z direction along field lines. At a given position the Alfven perturbations 
are torsional oscillations (i.e. oscillations are in the ^-direction, perpendicular to the background field). (B) 
Alfven waves propagating along a constant magnetic surface. The Alfven perturbations are within the mag- 
netic surface (yz plane) at the discontinuity, perpendicular to the background field (y direction), whereas the 
waves themselves propagate along the field lines (z direction). Density enhancements (e.g. in the form of 
spicules as observed by De Pontieu et al. 2007 or within prominences as seen by Okamoto et al. 2007) are 
visualised here as a yellow-red thin blob that follows the field lines. Vertical arrows indicate the magnetic 
field gradient increasing from left to right (reproduced from Erdelyi & Fedun 2 007) . 



exceed the critical twist for the MHD kink instability (e.g. Raadu| ( |1972| >). The z-axis of 
the cylindrical coordinate system (r, ip, z) is along the tube. The radius of the tube (ro) is 
constant. The unperturbed state and geometry of the model of a magnetically twisted tube 
is sketched in Figure [TO] All dependent variables inside the flux tube have index i, while 
quantities outside the tube are denoted with index e. In cylindrical equilibrium the magnetic 
field and plasma pressure satisfy the equilibrium condition in the radial direction: 

~r~ I Poi H 7. H = 0. (2) 

dr ^ 2^0 J Hor 

Here, the second term in the brackets represents magnetic pressure and the third term of 
Equation Q derives from magnetic tension due to the azimuthal component of the equi- 
librium magnetic field, Bo lip . Gravity is not included in the present analysis. The particular 
case of a uniformly twisted equilibrium magnetic field, Bo(r), of the form: 

B o« = (mn r R B ° iz) ' T°' ® 

1 (0,0,B 0ez) , r > r , 



is considered, where 



A 2 r 2 



1/2 



B Qiz =Bll- 2-^- I (4) 

in order to satisfy the pressure equilibrium governed by Equation Here A and B are 
arbitrary constants. This choice of magnetic field may represent solar atmospheric flux tubes 
with observed weakly twisted field components (see e.g. |Klimchuk et al.|2000[ >. 
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Fig. 10 The equilibrium configuration. The uniformly twisted magnetic flux tube is in an ambient straight 
and uniform magnetic field. At the boundary of the flux tube there is a jump in magnetic twist (reproduced 
from|Erdelyi & Fedun 2010). 



In a cylindrical coordinate system with the equilibrium state inhomogeneous in the radial 
direction only, £ = (£ r , £ z ) is the Lagrangian displacement vector denoting perturbations 
from the equilibrium position; the unperturbed mass density po and pressure po are constants 
(though not equal to zero) inside the flux tube, respectively; Bo = (0, Bo v (r), Bo z (r)) 
is the equilibrium magnetic field; no is the magnetic permeability of free space, p is the 
perturbation of kinetic plasma pressure and b is the perturbation of the magnetic field. Let 
us introduce the Eulerian perturbation of the total pressure: Pt := p + Brjb/^o- Because 
the equilibrium variables depend on radius only, for normal-mode analysis in cylindrical 
geometry we Fourier decompose the perturbed variables with respect to (p and z through, 

£, Pt ~ exp[i(fc,z + mip — uit)] . 

Here uj is the mode frequency, m is the azimuthal order of the mode, k is the longitudinal 
(axial) wavenumber. Then, after some algebra using the ideal linear MHD equations, one 
can find that £ r and Pt satisfy a system of two first-order ordinary differential equations 
(see e.g. |Appert et al.|1974||Sakurai et al.|1991^ , 

D ^H*) = - ^Pt ,and (5) 

dr 

D-^Pt = -C* 3 r£r - CiP T . (6) 
dr r 

The choice of variables above does not give much insight into the governing equation 
of Alfven waves. We may need to work on further for a more suitable choice of variables, 
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in order to easier identify perturbations that may be characterised as Alfven waves. We 
will do this in due course. Here, the coefficient functions D, C\, C 2 and C3 depend on 
the equilibrium variables and on the frequency w. In general, these coefficient functions are 
given as follows: 

D = Pq(oj 2 — ui\)C/t, (7) 
Ci = 2B^ / 4 B rn \ (8) 
Hor \ r J 

C 2 = u, 4 - (k 2 + ^pj C 4 , (9) 
V P0P0 dr V r I I \ p r 



45, 



2 

PvCa^ujI, (10) 



where 



and 



C 4 = (C§ + Vt) (lo 2 - u 2 c ) , (11) 



2 PO 2 2 2 2 

Cs = 7 — > Va = > B = B 0ip + B 0z , 

Po P0P0 

t m DilD 2 /b 2 C| 2 

r popo C| + 

Here Cs is the sound speed, Va is the Alfven speed, uia is the Alfven frequency and ujq is 
the cusp frequency. To obtain the remaining perturbed variables, i.e. £ z , b r , b v , b z and p 
in terms of £ r and Pt we will have 



2f 

P0^ ?r 


d n , 2Bo v b v . / B fe r 
dr /i r /io 


(12) 




= 1— P T - 1 - j (VBovO , 

r po Po r dr 


(13) 


2^ 
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sip i/s&z &r d 

= lfcP T 3-^Oz, 

Po Po dr 


(14) 






(15) 


b v 


= ifc (B z£<p — B 0v £ z ) — — (Boip£r) , 


(16) 


6z 


= — ' — (Boz£v> — B 0ip ^) j- (rBoz^r) , 

r r dr 


(17) 
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/id,., im ., \ 

= -IPO -T-( r Cr) H C V +*fc£z • 

\ r dr r / 


(18) 



From the above equations, the dominant Alfven perturbations are described by Equations l |13| l 
& {16}. Note, that the perturbations are coupled to the other wave modes, and one may re- 
cover the torsional Alfven wave only if the conditions m = and Bo v = (i.e. no equilib- 
rium twist) are concurrently satisfied. Straightforward algebra yields the following expres- 
sions for the perturbed quantities as functions of the radial component of the displacement 
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vector (£ r ) and the Eulerian perturbation of total pressure (Pt): 
D \ r p po J 

i/Bgpy ( „ 2 f 2 m /b r> 
977 04 Po^ i>0« 

poPo^D \ r \ r mo 

£z = -=: fcC 4 W P T + 

i/bBoz / '„ ,2 m / B n \ t 
popooj z D \ r r fio I 



k 



Po(^ 



rp (ui 2 - u)\) 



3bPt & 

Por 



2 \ 1 SbPt ^ £. 



Po r 



ld_ 
r dr 



(r-BozCr) , 



(19) 
(20) 

(21) 

(22) 

(23) 
(24) 



where 



mB 0z 
g B = kB { 



Oip ■ 



(25) 



Note, that the radial perturbed component of the magnetic field b r is determined directly by 
Equation (JT5J, In many solar applications it is more practical to use perturbations parallel 
(||) and perpendicular (_L) to magnetic field lines within constant magnetic surfaces. These 
perturbations will represent longitudinal and Alfven perturbations. For example, in the case 
of no magnetic twist (i.e. A = 0) and for axially symmetric perturbations only (i.e. m = 
0) the perpendicular perturbations of the magnetic field and velocity within the constant 
magnetic surfaces would describe the torsional Alfvem waves. In the geometry of the present 
paper such surfaces are concentric cylinders centred about the z-axis. We introduce the 
components of the displacement vector and the perturbed magnetic field parallel (£y, by) 
and perpendicular (£x> ^_l) to me background equilibrium magnetic field lines within the 
constant magnetic surfaces given by, 



(£y,by) = ({^,b v )B 0ip + (^,b z )B Qz )/B , 
(f_L,&x) = ((S<p,b v )BQ Z -(&,b z )Bo<p)/Bo. 



(26) 
(27) 



These components are related to £ r and Pt as, 
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Fig. 11 The dimensionless phase speed (V p h) as function of the dimensionless wavenumber (fcro) under (a) 
typical photospheric (i.e Vai > Cs c > Csi > Va c ) and (b) coronal conditions (i.e Vac > Vai > Csi > 
Cg e ) for an untwisted magnetic flux tube are shown at the (a) left and (b) right panels, respectively. Here 
(a) C Se = 0.75V A i, Vac = 0.25Vai, C* Si = 0.5V A i for the photospheric case and (b) C Sc = 0.25V A i, 
Vac = 2.5VAi. Csi = 0-5VAi for the solar corona. In the left panel three modes of the fast surface 
waves that are bounded within [Cgi, Cg c ] and the infinite number of the slow body waves bounded within 
[Cxi, Csi] are plotted. The slow kink, sausage, etc. surface waves are very close to each other (just under 
Cxi). In the right panel sausage (m=0), kink (m=l) and fluting (m=2) modes are shown. The infinite number 
of the fast body waves are bounded within [Vaj, Va c ] and the infinite number of the slow body waves are 
bounded within [Cxi, Csi]. m = ij, where the i refers to the mode (sausage, kink, etc.) and j refers to the 
j-th branch of the zeroes of eigenfunctions in the radial direction. Only three branches of each mode of the 
infinitely many slow and fast body waves are plotted. 



These governing equations, in various alternative formats, were derived previously in 
various contexts by e.g. Hain & Lust ( 1958), Goedbloed ( 1971 ), Appert et al. ( 1974), Goossens 
< [!99lT ), |Sakurai et al.| ( fl99"l| l and |Erdelyi & Fedun| i |2006||2007||2010] >~ 



5.3 Derivation of the General Dispersion Equation 

Applying the boundary conditions to the inside and outside solutions of the governing equa- 
tions yields the required general dispersion relation that will determine the existing wave 
modes of the system. The continuity condition for the radial component of the displace- 
ment, 

Cir| r=ro = £er| r=ro , 
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and the continuity of pressure, 



B, 



Oitp 



fofo 



6, 



Te\r=r 



provide the desired dispersion relation: 



r K m {m 0e ro) 
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/ 1 u M(a,o,x ) 



(29) 



The dash denotes the derivative of the Kummer's function evaluated at x = xrj, where 

n cil/2, 2 2 ,, m 

x = -E 1 k a r Q . (30) 

Equation \29\ is the general dispersion equation for linear MHD waves in a magnet- 
ically uniformly twisted compressible flux tube embedded in a compressible and uniform 
straight magnetic environment. To the best of our knowledge, it is not possible to solve ana- 
lytically Equation ( |29| l in general terms. However, at least certain checks can be carried out, 
e.g. applying the incompressible limit should recover the dispersion relation of Ben nett et| 
1999| > or Erdelyi & Fedun ( 2006 1. Another limiting case is Bq v — > where the well 



al. 



known dispersion relation Equations (8a, b) of jEdwin & Roberts| ( |1983| l must be recovered. 
These particular cases of analytical checks reassure us that the derivations are consistent. 
Although we would prefer further analytical insight, such approach seems to be limited by 
the lack of existing suitable asymptotic expansions for Kummer functions M (a, b, xo) and 
M (a, b, xo). We shall here only recall a typical graphical solution of the dispersion relation, 
revealing all the linear MHD modes. 

In Figure [TT] the results of calculations of the dispersion curves derived from Equation 
\29\ are presented. One can observe here the very rich structure of linear MHD modes. 
Note, the Alfven waves are just a single line, as all the modes have the same phase speed. 
The modes can often only be distinguished from other magneto-acoustic (e.g. kink) modes, 
if one carefully determines the polarisation. This can be rather hard, given the currently 
available observational tools. 



6 Concluding remarks 

There are at least three fundamental questions that must be resolved in order to understand 
the heating of the solar (and stellar) atmosphere: Where is the energy generated? How does 
the generated energy propagate from the energy reservoir to the solar atmosphere, spanning 
from the warm chromosphere to the corona? How does the transported energy efficiently 
dissipate in the solar atmosphere in order to maintain its multimillion-Kelvin temperature? 
There are numerous examples of slow and fast MHD waves observed in the solar atmo- 
sphere. However, their energy density does not seem sufficient for coronal heating or pow- 
ering the solar wind. On the other hand, among MHD wave theorists, Alfven waves still 
remain the "Holy Grail" of coronal heating, solar wind acceleration and solar magneto- 
seismology (Banerjee et al. 2007jl. 

The applicability of MHD wave research to magneto-seismology is a relatively recent 
development. The main idea here is to determine information about the MHD waveguide 
in which the waves propagate. This is analogous to the techniques of seismology, a popular 
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method of investigating the sub-surface structure of the Earth. The difference here is that in 
a magnetised plasma there is a plentiful abundance of waves in existence, offering us various 
aspects of seismology. The frequency (or period) measurements of MHD waves are often 
used to determine density and magnetic stratification, both in the longitudinal and transver- 
sal directions of the waveguide. More recently, however, studies of the spatial distribution 
of waves and their intrinsic characteristics have been suggested, offering an independent 
method to determine diagnostic information about the waveg uide ( |Erdelyi & Verth|2007| 
|Verth et "aL]|2007| l. Magneto-seismology has a novel application to the reconstruction of 
MHD waveguides. This is perhaps an even more important aspect compared to the others. 
The idea here is that by studying the properties of MHD waves and oscillations, their dis- 
tribution in the waveguide (assuming here an adequate spatial and temporal resolution) may 
help construct detailed magnetic surface topology. The idea has been outlined in a recent 
paper by |Verth et aT1pjT0] i and |Fedun et~aT1, l |201 lb) . 

There have been some suggestions that transversal waves may be considered Alfvenic 
in nature providing that they are incompressible, exhibit no intensity fluctuations and dis- 
play periodic displacements with the magnetic tension as the restoring force (Goossens et 
|al.|20 09). The terminology of Alfvenic is, however, not entirely clear as such waves do not 
normally exist as an eigenvalue of the linear MHD equations in structured plasmas: on the 
other hand many magnetic oscillating structures in the corona and below can not be de- 
scribed by pure linear MHD eigenmodes. Given the ambiguity, it may be better to refer to 
these waves as transversal ( [Jess et al. 2012 1. Further, refined observations, with the use of 
e.g. spectro-polarimetric measurements, may be able to clarify the wave mode identity with 
greater accuracy. 

According to Cirtai rTet al.| ( [2007 \ and |De Pontieu et al.| < |2007| ), signatures of Alfven 
waves were observed with the use of the XRT and SOT instruments, respectively. Because 
these jets are increases in local density, if their length is short relative to the wavelength of 
Alfven perturbations, they can slide along magnetic field lines, and act as test particles (see 
the yellow-red blob in Fig. [9}}) in the plasma that has almost no resistivity, just like bobsleds 
would run and slide along a prebuilt track. However, these observations also raise concerns 
about the applicability of the classical concept of a magnetic flux tube in the very dynamic 
solar atmosphere, where these sliding jets were captured. In a classical magnetic flux tube, 
propagating Alfven waves along the tube would cause torsional oscillations (Fig. [9^). In this 
scenario, the only observational signature of Alfven waves would be spectral line broaden- 
ing, as observed by Jess et al. ( 20 09|l. On the other ha nd, if these classical flux tubes did 
indeed exist, then the observations of De Pontieu et al. (2007) would be interpreted as kink 
waves (i.e. waves that displace the axis of symmetry of the flux tube like an S-shape). More 
detailed observations are needed, so that a full three-dimensional picture of wave propaga- 
tion could emerge. 

Numerical simulations of a two-dimensional stratified VAL-atmosphere model (Ver- 
|nazza et al. 19811, embedded in a horizontal magnetic field to mimic solar prominences and 
driven by global photospheric motions, will also result into observational signatures very 
similar to those reported by Okamoto et al.| ( |2007| l, and in Figure 2 and the online movie 
SI of Er delyi & Fedun| l |2007 1. In particular, note the similarity between the time-distance 
plot (Fig. 2a) of the illustrative forward modelling simulations in Er delyi & Fedun| ([2007 1 
and Figure 3 of Okamoto et al. ( 2007 1. However, in the case of this numerical forward mod- 
elling experiment the oscillations found are clearly not Alfvenic, but are rather transversal 
kink oscillations. 

It is worth mentioning what may be some key directions where MHD, in particular, 
Alfven wave research may develop. Among the many possibilities is the development of 
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Alfven (and in parallel other MHD) wave theory in dynamic plasma. Observations indicate 
that waveguides change, often considerably, during the time that the oscillations reported are 
guided by them (Aschwanden & Terradas 20081. As long as the characteristic time scale of 
the changes in the background plasma are small or comparable to the time scales of waves 
crossing the structure, some initial theoretical progress has already been made (Mor ton &| 
|Erdelyi|2009l[20T0tporton et al.|201Ql[20TTt|Ruderman|201 la|b) l. It is a real challenge to 
develop wave theory in dynamic plasmas. Another fascinating aspect of developments will 
be the advancement of (non)linear Alfven and magneto-acoustic waves where the geometry 
and structuring of the waveguide is more complex but more realistic. Significant techno- 
logical improvements associated in High Performance Computing will allow detailed MHD 
wave studies of the coupled solar interior-atmosphere system as a whole. Such studies will 
shed ample light into the true nature of Alfven and other MHD wave dynamics and heat- 
ingjBogdan e t al.|2003H"Carisson & Bogdan|2006||Erdelyi et al.|2007[|Fedun et al.|2009l 
Khomenk oet al.|2008[|Fedun et al.|201 la^ . 

Last, but by no means least, we will provide some observational goals that will hopefully 
assist the detection, disambiguation, and study of solar oscillations. High-cadence imagers, 
such as the ROSA ( [Jess et akp OlO) instrument, will need to be employed with high spatial, 
temporal, and spectral resolution facilities, such as the IBIS (Cavallini 2006 ) and CRISP 
(Scharmer et al. 2008) two-dimensional spectro-polarimeters. Simultaneous observations 
between the ground-based instruments and state-of-the-art space-borne facilities such as 
SDO, and the planned IRIS (Interface Region Imaging Spectrograph) can provide the nec- 
essary intensity, magnetic-field, and Doppler signals required to unveil the specific wave 
modes as they travel through the solar atmosphere. The Helioseismic and Magnetic Imager 
(HMI) on SDO provides dopplegrams, continuum filtergrams, line-of-sight and vector mag- 
netograms at the solar surface, thus allowing the sites of photospheric oscillations to be 
identified (albeit at a lower spatial resolution). The Atmospheric Imaging Assembly (AIA) 
provides the coronal magnetic field configuration, which can be compared to that obtained 
by HMI. IRIS in particular will obtain high spatial (250 km) and temporal resolution (1 sec) 
spectra and images in ultraviolet wavelengths, focusing on the upper chromosphere and 
transition region. The combination of simultaneous imaging and spectroscopy provided by 
IRIS will allow us to disentangle the velocity and intensity signal of torsional, kink and 
sausage wave modes as they travel through the complex interface between the photosphere 
and the corona. High resolution limb spectro-polarimetry of chromospheric structures such 
as spicules and prominences remains a major challenge. The low photon statistics at the 
solar limb, limit the use of Adaptive Optics and the number of photons required for the mea- 
surement of accurate Stokes profiles. These difficulties will be overcome with the upcoming 
Advanced Technology Solar Telescope ( Keil et al.|2010| >, with a 4 m primary mirror, will 
allow structures in the Sun's atmosphere to be examined on unprecedented spatial scales of 
r;20 km. 
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